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Recent progress in experimental studies of low-dimensional systems with strong spin-orbit coupling 
poses a question on the effect of this coupling on the energy spectrum of electrons in semiconductor 
nanostructures. It is shown in the paper that this effect is profound in the strong coupling limit. 
In circular quantum dots a soft mode develops, in strongly elongated dots electron spin becomes 
protected from the effects of the environment, and the lower branch of the energy spectrum of 
quantum wires becomes nearly flat in a wide region of the momentum space. 
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I. INTRODUCTION 

Quantum dots are basic elements of the systems aimed 
in semiconductor based quantum computation based 
on qubits encoded in electron spins^— Electrical op- 
eration of such qubits is possible through the Heisen- 
berg (exchange) mechanism, 4-6 or spin-orbit coupling, 
or both. 7 Effective spin-orbit coupling in quantum dots 
originates from two different mechanisms. First mech- 
anism is based on the celebrated Thomas term cou- 
pling electron momentum p to Pauli matrices er. In 
two-dimensional (2D) systems it is usually represented 
through two invariants linear in the quasimomentum k, 
known as the Rashba and Dresselhaus terms, reflecting 
the structure, bulk, and interface induced asymmetry. 9 
Second mechanism couples <x to the electron coordinate r 
through the inhomogeneous external, exchange, or hyper- 
fine magnetic field . 10 ' 11 Electric dipole spin resonances 
(EDSR)2ii2~— driven by these mechanisms have been 
observed experimental!}*^— in different semiconductor 
nanostructures. Electric fields generated by precessing 
spins are an inverse effect but nediated by the same 
mechanisms i^i 

Effect of spin-orbit coupling on the spin dynamics in- 
side a dot can be quantified by the ratio of the dot size d 
to the spin-orbit length £so- For GaAs, spin-orbit cou- 
pling is relatively weak, with £so ~ 35 /im, hence, with 
d ~ 100 nm, d/£ S o ~ 10 -3 < 1 J£ Currently, interest is 
shifting in the direction of semiconductor dots from the 
materials with much stronger spin-orbit coupling such as 
InAs and InSb; for them lso ~ 300 nmJ£ For all materi- 
als with d/£so 1) spin-orbit coupling influences mostly 
the spin dynamics and has only a minor effect on the 
shape of the electron cloud inside the dot. 

However, recently new classes of systems emerged, 
based on the chemical elements with high atomic num- 
bers, where spin-orbit coupling is even much stronger, up 
to the atomic scale. They are considered as new potential 
materials for spintronics,— in particular, for the Datta- 
Das spin transistor^ If to estimate spin-orbit energy 
in the Hamiltonian as asok, in such systems aso 1 
eV A. These are noncentrosymmetric bulk BiTel with 
a so = 3.85 eV A,2i Te-terminated surface of BiTeCl 



with a = 1.78 eV A and the minimum of the surface 
electron band well inside the bulk gap, 25 HgTe quan- 
tum wells with fc-dependent a reaching the value a < 1 
eV Aj2£ and interfaces between insulating perovskite ox- 
ides LaA10 3 /SrTi0 3 with a so <5x 10~ 2 eV k^r^. 
to name only a few (topological materials are outside 
the scope of the paper). For nanocrystals of high aso 
compounds, and quantum dots in them, an opposite in- 
equality, djlso ^ 1, should be fulfilled, with the charge 
distribution strongly influenced by spin-orbit coupling; 
for quantum dots (or nanocrystals) of BiTel and BiTeCl, 
one expects d/£so ^ 1- The focus of the present pa- 
per is on the electronic structure of strongly spin-orbit 
coupled 2D nanostructures, but still in the range of the 
applicability of the envelope function approximation. 



II. PARABOLIC ALLY CONFINED DOTS 

A model Hamiltonian for a lateral quantum dot with 
parabolic confinement is 
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It is convenient to introduce dimensionless units using 
%ujq as a unit of energy, E = frcooe, and £q = V"7 muj o as 
a unit of length. Then, with r = lop, k = k/£q, £so — 
Ti 2 jmasoi an d a = ^o/^so: the Hamiltonian reduces to 



h = 



a{cr x k) z 



El 

2 



(2) 



The energy spectrum of a free particle [last term in 
Eq. @ omitted], consists of two branches eo( K ) = k 2 /2 — 
an with a minimum e a = —a 2 /2 achieved at a circle 
of the radius of K a = a in the momentum space. For 
finding the energy spectrum of a quantum dot near its 
bottom in the strong spin-orbit coupling limit, a ^> 1, 
it is convenient to solve the eigenvalue problem for a 

u 



two-component spinor tp 
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in the momentum 



u + ian-V = 0, 
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v — ian+u = 0, (3) 



with k± = «e ±l >. Therefore, u{p,<p) = u(p)e^' -1 / 2 ^, 
= -y(p)e i (J+ 1 / 2 )^, with half-integer j > 1/2; levels 
with j < possess the same energy due to the Kramers 
symmetry. Equations for the radial parts (u(p),v(p)) are 
simplified when rewritten for their linear combinations 
w± = uztiv 
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It is seen from Eq. ([3|) that near the spectrum bottom 
| itf— | 3> |i0+|i an d w+ w (j/2o! 4 )u>_. Because in what 
follows the expansion will go in the parameter j/a <C 
1, Eqs. (@]) decouple. In the same spectrum region, the 
equation for w- reduces to 
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First term, enclosed in parentheses, is a harmonic os- 
cillator related to the radial motion with the spectrum 
(n = (n+ 1/2) — a 2 /2, n > 0. Second term produces, in 
2nd order of the perturbation theory, an n-independent 
level shift Ae = — l/8a 2 ; it plays no significant role. 
Third term describes a soft mode related to the angular 
degree if freedom with a typical level spacing of 2] /or; 
first gap, separating the ground and first excited states, 
equals 2/ a 2 <C 1. Finally 
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with E a = —Tiuioa 2 /2, n > 0, and j > 1/2. Small sep- 
aration between the energy levels with the same n but 
different values of j hampers trapping of a single elec- 
tron into a well defined quantum state. 

For comparison, in absence of spin-orbit coupling, a — 
0, the Hamiltonian of Eq. Q describes a two-dimensional 
harmonic oscillator. Its ground state is well protected by 
gaps of which the smallest one equals to Tiujq. 

The above results prove that, in axially symmetric sys- 
tems with a strong spin-orbit coupling, a» 1, confine- 
ment in quantum dots is hindered by the motion along 
the angular degree of freedom, and achieving a well iso- 
lated ground state requires freezing angular dynamics 
by breaking the axial symmetry. Angular dynamics is 
completely suppressed in a one-dimensional setup with a 
Hamiltonian 
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that follows from Eq. ([2} by suppressing first Cartesian 
coordinate in p = {£,rj). It possesses the harmonic os- 
cillator energy spectrum e„ = (n + 1/2) with the ground 
state eigenfunctions 



ipi(r]) 
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forming a Kramers doublet; (cri,(T2) are two first Pauli 
matrices. Hence, the energy spectrum is well gapped 
(with the gap AE = huio), the electron density is a 
smooth function of the coordinate and is trapped near 
the center of the parabolic well at the scale Ay ~ £o, 
while the mutual phase of the spinor components oscil- 
lates with a frequency controlled by the spin-orbit cou- 
pling constant a. Fourier images of (ipi,ip2) are peaked 
at k v = ±ia. 

However, as is explained in Sec. IIVI below, the Hamil- 
tonian of Eq. ([7]) is not a proper model Hamiltonian for 
a quasi-one-dimensional quantum dot in a strongly spin- 
orbit coupled material described by the Hamiltonian of 
Eq. @. In particular, the Fourier images of wave func- 
tions are peaked at k v = ±a, hence, they imply a state 
that is a double dot in the momentum space. The origin 
of this striking difference comes from the fact that two 
first terms of the Hamiltonian h of Eq. ([7]) with a = const 
do not describe the dispersion of the lower branches of the 
energy spectrum of a quantum wire that are flat in a wide 
region of the momentum space \k v \ < a, see Fig. 1. Un- 
usual properties of quasi-one-dimensional quantum dots 
in strongly spin-orbit coupled materials are discussed at 
the end of Sec. M\ 



III. HARD- WALL CONFINED DOTS 

Hard-wall confinement provides an independent insight 
onto the same problem. It has been already investi- 
gated by Bulgakov and Sadreev^i and Tsitsishvili et ali^ 
Therefore, technical details will be omitted in what fol- 
lows, and we focus on the large a so limit; the Hamil- 
tonian is the same as in Sec. IIVI below. After choosing 
the radius R of the confinement region as the unit of the 
length, introducing dimensionless quantities p — r/R, 
e = E j (h 2 1 mR 2 ) , a = R/£so, an d separating the an- 
gular variable <j) similarly to Sec. [Til equations for the 
eigenspinor components (u, v) are 

[p 2 d 2 + pd p - (j-l/2) 2 + 2ep 2 ]u 
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[p 2 d 2 + pd p - (j + l/2) 2 + 2ep 2 

i-1/2. 



+ 2ap 2 [d p - J 1 ]u=0. 

Due to the recurrence relations^ 

[d p ± (j ± 1/2)/ p] Jj±i/ 2 (M = ±kJ jt1/2 (k P ), 
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solutions of Eqs. ([9]) are Bessel functions u(p) = 
aJj_i/2(Kp), v(p) = bJj +1 / 2 (np), with half-integer val- 
ues of the angular momentum j > 1/2. The hard- wall 
boundary condition u(l) = v(l) = results in the fol- 
lowing equatio n 31 ' 32 for the energy eigenvalues e 

Jj-l/2(K+)J j+1 /2(K-) = J j - 1/2 {K-)J j+1 / 2 (K + ), (11) 

with 



j. The ground state is protected by gaps of which the 
smallest one is approximately equal to 4A5h 2 /mR 2 . 

Equation (fT5|) is valid only near the spectrum bottom, 
£nj - e Q <C e Q . Numerical analysis of Eq. (|11[) shows 
that with increasing j its roots move up, and the total 
number of negative roots decreases. Finally, at j ~ a, 
there remains a single negative root < that with 
increasing j asymptotically approaches zero, e\ t j 0. 



= a ± yc 



2e. 



(12) 



Equation ([TT]) was derive d 31 ' 32 under the condition of the 
existence of two spectrum branches, therefore, its triv- 
ial solution k + = k_, corresponding to the free-particle 
spectrum bottom, should be omitted. 

Below the conical (Dirac) point of the spectrum, for 
—a 2 /2 < e < 0, both k± > 0. Near the spectrum bot- 
tom, k± w q > 1, energy spectrum can be found from 
the k> 1 asymptotic expansion of Bessel functions^ 3 . 
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In the leading order in 1/k <C 1, with e close to the 
spectrum bottom e a = —a 2 /2, when e — e a <C a and 
j <C a, Eq. reduces to 



jy/2(e - e a ) { cos(2a - jn) + j cos [2 ^2(e - e Q )]} 

= a 2 sin[2V2(e-e Q )]. (14) 

In the leading order in 1/a <C 1, solutions of Eq. (fT4|) 
are 2-^/2(e — e a ) ~ 7rn, with n > 1. In the next order, 
for y/Jn <C a, the energy spectrum near its bottom E a — 
—H /(2m£ 2 ;o ), expressed in dimensional units, is 



E„,j — E a 



h 2 n 2 n 2 
8mR 2 
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with integer n > 1 and half- integer j > 1/2. It is seen 
that there are two quantum numbers. 

Integers n quantize the radial motion similarly to a 
particle confined in a hard-wall box in absence of spin- 
orbit coupling. This coupling does not change the scale 
of the level spacings of the radial motion that is set by 
the radius R of the confinement region. It is seen from 
the form Jj±\/2(np) of the solutions of Eq. © that eigen- 
functions are oscillating functions of p. 

Half-integers j describe angular motion. Remarkably, 
second term in the braces including the angular quantum 
number j shows oscillatory dependence on a. Similarly to 
Eq. ([6]), the angular mode is soft with the level separation 
small in 1/a 2 . 

For comparison, in absence of spin-orbit coupling, a = 
0, the energy spectrum is described by a Bessel equa- 
tion similar to Eq. ^ but with integer angular momenta 



IV. QUANTUM WIRES 

Interest in the energy spectrum of single- and multi- 
channel spin-orbit coupled quantum wires increased re- 
cently after they have been proposed as a platform for 
solid state realization of Majorana fermions^~— 

For a hard- wall confined quantum wire of a width Y, 
—Y/2 < y < Y j 2, it is convenient to choose Y as a unit of 
length and Ti 2 jmY 2 as a unit of energy, E = (h 2 /mY 2 )e. 
Then, with dimensionless a defined as a = Yjlgoi the 
Hamiltonian is 



h = — + a(a x k) z . 



(16) 



In Eq. (|16|) . r and k are dimensionless coordinate and 
momentum, respectively. 

For energies below the conical (Dirac) point, e < 0, 
eigenfunctions and energy of a free particle are 



^ k(r) ~ ^2 ( Kk x + ik y )/k 1 ' 



;ik - r) ,e (fc) = Y~ ak - 
(17) 

It is convenient to define momenta 



k±(e) 



2c 



(18) 



at the external (internal) arcs of the energy spectrum. 
Then, for k x < fc_, eigenfunctions of confined states can 
be chosen as linear combinations of four functions ^k(r) 
with 



k±(e,k x ) = Jk 2 ± (e)-kl 



(19) 



and —k^(e,k x ), obeying the boundary conditions 
ipk(x, ±1/2) = 0. After some algebra, the eigenvalue 
problem reduces to 

(k 2 + 2e) smk~ sinfc^ + fc~fc+(l — cos k~ cos fc+) = 0. 

(20) 

Equation (|2T)|) can be solved exactly for k x = 0. Indeed, 
in this case kyky = 2|e| and, for e < 0, Eq. ([20]) reduces 
to cos(fc + — fc_) = 1. Its solutions are 



fc+(e) — fe-(e) = 27m, n > 1. 



(21) 



The value n = is excluded because it would result in 
fe+(e) = fc_(e) and decreasing the number of linearly in- 
dependent basis functions. Squaring Eq. (f2"Tj) leads to 



e n (k x = 0) 



n > 1. 



(22) 
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Therefore, for k x — the energy spectrum reduces to 
the spectrum of the textbook quantum box problem, and 
spin-orbit coupling manifests itself only through lowering 
the spectrum bottom to e a = —a 2 /2. Energy levels of 
Eq. (1221) are Kramers degenerate and split with increasing 
k x as seen in Fig. 1. 

The spectrum described by Eq. (|2T))) extends until 
k x = fc-(e), which corresponds to e = k x /2 — ak x . This 
boundary is shown by a heavy line in Fig. 1 and describes 
a free particle moving in the x direction. To the right 
from this line, propagating waves exp(±ifc~y) change to 

evanescent waves exp(±K y y) with n y = k 2 — Ari(e), 
and Eq. (f^TJj) should be substituted with 

(k x + 2e) sinh n y sin ky + n y ky ( 1 — cosh n y cos ky) = . 

(23) 

Energy spectrum for a = 10 is shown in Fig. 1. At 
k x = 0, there are three Kramers degenerate energy levels 
with energies e < that split with increasing momentum 
k x . For the n = 1 levels, this splitting is very small until 
k x ss 8, and there are two level intersections, near k x = 
3.0 and k x = 7.4. These are real intersections rather than 
avoided crossings because the Hamiltonian of Eq. (Tl7j]) 
possesses a xz symmetry plane and two n = 1 branches 
have opposite parity with respect to the reflections in this 
plane. The splitting only increases for a > 8 when the 
momentum approaches the value k x = a where the free- 
particle dispersion law eo(k x ) — k x /2 — ak x reaches its 
minimum e a = —a 2 /2. At k x = a, deviation of the lower 
branch from e Q is only e(k x = a) — e a m 105/2a 2 . The 
absolute minimum of the lower branch e m i n w e Q +0.411 is 
achieved slightly to the left from this point, at k x sa 9.442. 

The left arc of the free-particle dispersion curve eo(k x ) 
cuts the plane into two parts. To the left from it quantum 
states are described by solutions of Eq. (|2T)]) and to the 
right by solutions of Eq. (f23|). 

Flat shape of both n = 1 curves in Fig. 1, together 
with the position of the spectrum minimum near k x = 
a ^> 1, suggests unusual properties of the quantum dots 
prepared from large-a quantum wires. We discuss them 
below at a qualitative level. 

The curvature of the lower n — 1 branch near its min- 
imum is close to that of the eo(k x ) curve, i.e., close to 
unity. For localizing a particle close to the minimum, the 



width of the quantum state in the k x space should be 
Ak x <C 7r as it follows from Eq. (|2"2"1) . Hence, according 
to the uncertainty principle, the length of the dot X (in 
dimensional units) should satisfy the condition 

X > Y/ir. (24) 

If this condition is satisfied with sufficient accuracy, the 
ground state of the dot possesses highly unusual proper- 
ties. This is a Kramers doublet comprising the states 
centered at k x — ±q with spins polarized parallel to 
y in opposite directions; it resembles a double dot but 
in the momentum rather than in the real space. Tran- 
sitions between these states are only possible through 
tunneling in the k x space under a barrier of the hight 
of nearly 7r 2 /2 and the width of about 2a. Therefore, 
we encounter a unique situation when strong spin-orbit 
coupling protects electron spin against decoherence. Ap- 
plying a magnetic field B || y splits the ground state into 
a Zeeman doublet. Remarkably, the intensities of both 
the electron spin resonance and EDSR transitions be- 
tween the components of the doublet should be strongly 
suppressed because the tunneling process is involved. 39 
Therefore, in this unique regime electron spin becomes 
well protected against low-frequency electric and mag- 
netic noise by strong spin-orbit coupling. 

With decreasing X, wave function spreads fast in the 
k x space, mostly in the direction of small k x values, in- 
volving the upper n — 1 branch, and than also next 
branches of the energy spectrum. As X approaches Y, 
the shape of the dot becomes close to circular, and a soft 
mode discussed in Sees . HT1 and IIIII becomes pronounced. 

In conclusion, with the increasing strength of spin- 
orbit coupling, when the spin-orbit length becomes 
shorter than the nanostructure size, this coupling influ- 
ences heavily the energy spectrum and electron density 
distribution. In the strong spin-orbit coupling limit, cir- 
cular quantum dots acquire a soft mode, dispersion of the 
lower mode of quantum wires becomes flat in a wide re- 
gion of the momentum space, and electron spin becomes 
protected from the effects of the environment in strongly 
elongated quantum dots. 
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FIG. 1. (Color online) Energy spectrum of a strongly spin- 
orbit coupled quantum wire. Energy e(k x ) plotted vs mo- 
mentum k x for the coupling constant value a = 10; dimen- 
sionless units. Only negative part of the spectrum, below 
the conical point, is shown. Energy levels at k x — are 
numerated by a quantum number n > 1. Two branches 
emerging from each k x = quantum level are shown by 
full and dashed lines. Heavy (red) line shows free-particle 
energy eo(k x ) — k x /2 — ak x . Notice week dispersion of both 
n = 1 branches. 



